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Abstract 


Manipulators, 

7J * — 


A new class of asymptotically stable adaptive control laws is introduced for application to 
the robotic manipulator. Unlike moat applications of adaptive control theory to robotic 
manipulators, this analysis addresses the nonlinear dynamics directly without approximation, 
linearization, or ad-hoc assumptions, and utilizes a parameterization based on physical (time- 
invariant) quantities. This approach is made possible by using energy-like Lyapunov functions 
which retain the nonlinear character and structure of the dynamics, rather than staple 
quadratic forms which are ubiquitous to the adaptive control literature, and which have bound 
the theory tightly to linear systems with unknown parameters. It is a unique feature of 
these results that the adaptive forms arise by straightforward certainty equivalence adaptation 
of their nonadaptive counterparts found in the companion to this paper (i.e., by replacing 
unknown quantities by their estimates) and that this simple approach leads to asymptotically 
stable closed-Loop adaptive systems. Furthermore, it is emphasized that this approach doe3 not 
require convergence of the parameter estimates (i.e., via persistent excitation), invertlbility 
of the mass matrix estimate, or measurement of the Joint accelerations. 


1. Introduction 

In past years, many papers have appeared on the application of adaptive control theory to robotic 
manipulators (cf., [2]-[7], and Hsia [8] for overviewVT It is a general property of adaptive designs 
based on Lyapunov's Direct Method, that the Lyapunov /function is chosen as a simple quadratic type, well- 
known and well studied in the standard adaptive con/rol literature [12] [13]. However, this particular 
Lyapunov function was originally motivated for applications to the standard adaptive control problems 
(i.e., linear systems with unknown parameters), and not for nonlinear dynamical systems. Hence, applications 
of standard adaptive control techniques to robotic manipulators invariably require rhe dynamics to be considered 
as linear. This in turn, requires the use of ad-hoc assumptions and/or analysis techniques including 
1) treatment of position dependent quantities fis unknown constants, for which they must be assumed to vary 
slowly with time; 2) linearization of the system about some local operating point-valid only for small 
excursions from nominal; 3) the use of linear decoupled models for the links, which neglects nonlinearities 
and crosscoupling effects; and 4) neglecting the nonlinear and time-varying dynamics completely by assuming 
the plant la linear. Hence, stability results based on these assumptions are questionable, and a rigorous 
proof of stability for adaptive control of robotic manipulators remains unresolved. 

A recent exception to the above criticism is due to the work of Craig, Hsu and Sastry [9]. Here, a 
useful “linear in the parameters" formulation is exploited to simplify the analysis, and to demonstrate global 
convergence of an adaptive version of the computed-torque control law - without approximation to the nonlinear 
dynamics. However, the resulting adaptive controller requires the invertlbility of the mass matrix estimate 
(which is not guaranteed a-priori), and measurement of the Joint accelerations (which is generally unavailable). 

It is suggested In (9J, that the former can be handled by projecting parameter estimates into known regions of 
parameter space for which the mass matrix Inverse exists, and in which the true parameters are required to lie. 
However, knowledge and calculation of such regions is not straightforward and appears to be a weakness of the 
method. 

In this paper, the "linear in parameters" formulation of [9] is used in conjunction with a different 
Lvapunov function. Here, the choice of Lyapunov function is more closely related to the energy of the system, 
aud better retains the nonlinear structure and character of the dynamics. In addition, many problems associates 
with adapting the computed-torque control law directly are avoided by making use of the new class of exponentially 
stabilizing controllers introduced in [1]. Although these controllers are very similar in form to the computed 
torque method, they have many advantages in the nonadaptive case (cf., [1]), and have the unique property that 
they can be made adaptive by using a straightforward certainty equivalence approach (i.e., by replacing unknown 
quantities by their on-line estimates). Furthermore, the class of adaptive systems defined in this manner 
can be shown to be asymptotically stable without approximation to the nonlinear manipulator dynamics. This 
approach does not require convergence of parameter estimates (i.e., via persistent excitation), invertlbility 
of the mass matrix estimate, or measurement of joint accelerations. 

In the most recent literature (i.e., preprints, conference papers, etc.) there appears to be other work 
currently taking place which combines the linear in parameters formulation with a new Lyapunov function [10], 

[11]. Although this work is very new and is evolving very rapidly, we will try to contrast our results where 
possible, and provide an overall perspective. 
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The format of the paper la at followa. In Sac. 2 the raaulta of [1] ara reviewed and summarized aa 
required for treatment of tha adaptive control case In Sec. 3, asymptotic stability ia proved for the claaa 
of lyactni arising from certainty equivalence adaptation of tha control lava in [lj. Slightly tangential to 
the main thruat of the paper la tha analyala in Sec. 4 of the adaptive computad torque mathod. Since 
the coaputed- torque control lav ia widely establlehed In the literature, and widely applied in practlcm. 

It la uaeful to apply the tcchnlquaa developed herein to ace to what extent it can be aade adaptive and to what 
extent a tab ill ty can be guaranteed. In Sec. 5, aeveral renarka are aade pertinent to the new adaptive designs, 
and conclualona are given in Sec. 6. 

2. Background and Notation 
2.1 Manipulator Dynaaica 

The well-known Lagrange-Euler equations of notion for the n-joint aanlpulator la given as follows. 


M(q 1 )q 2 - - C(q lt q ? ) - k(q 1 ) + u 


l* M 2 y 


where 


C(q l’ q 2 ) " l(e i q 2 T V q l ))T ' 1 (, i q 2 T J 

^ unit vector 

3M(q i> a th 

Vi* “ 1^— - : *>ii “ 1 component of q^ 


k(q^) • gravity load 


( 2 . 1 ) 

( 2 . 2 ) 


(2.3) 


Here, ucR n la a generalized torque vector, q^, qj, q 2 cR n are generalized Joint position, velocity and 
acceleration vector, (e.g., ia an angle or a distance for a revolute or prismatic Joint, respectively, 
M(qi)cR nxn is the syraaetrlc positive definite mass inertia matrix; C(qj,q 2 )cR n la the Coriolis and centri- 
fugal force vector; and k(q^)cR n is the gravitational load vector. 


2.2 Some Useful Identities 


Let the following notations be detlned, 

W* 5 - I K^)* ej 

i"l n 

M( q i. q 2 ) - 57 M( qi ) \\ 

J(q 1 ,z) - l [(e 1 z T M 1 (q 1 ) - (e 1 z T M 1 <q 1 ) ) T ) 


N ' q i - q id * *2 ' q 2 - q 2d 


q ld* q 2d " desire<1 J OAnc position and velocities respectively (q^ ■ q^) 


T f l 


rCqj^.qj.q^) • dq 2 [j M(q 1> q 2 )Aq 2 - C(q 1> q 2 )q 2 ) 


Using the above notation, the following Identities are quoted from [i] without proof. In these 
identities, x, y and z are used to denote arbitrary vectors of appropriate dimension, 

Identity 1 

M ^l'^2^* " ^D^l* 2 ^ whjer « vector z la arbitrary 


Identity 2 

C(q r z)z - y (M D (q 1 .z) - J(q r z))z 
Identity 3 

Jtqj.z) * M D T (q 1 ,z) - M I) (q 1 ,z) 
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2.3 Important Leon* 

In this section, a useful l«m le reviewed, quoted directly without proof fro* fl]. For convenience, 
thle result will be alternstively referred to ee the g-Ball Tews due to the aethod used to prove it. 

Ltnu 2-1 (0-Bell Lease) 

4 

Given e dynamical systea 

*1 * f i* x r*“ x lT t * * ^t*" 1 • 1 i° 

Let fj's be locally Llpschltx with respect to X},...,x H uniformly In t on bounded intervals and continuous 
in t tor t > 0. Suppose a function V:R n l x **' xn N x * *+ is given such that 

r T 

V(x 1 ,..., V t) - ^ l x t P 1J (x 1 X J ,t)x J , 

V is positive definite in % 1 x^ 

*<«! V° - ' I (o i - I > lj ll* i (‘)|| k «)||x 1 (t)|| 2 (2-*) 

UI 1 J el 2 t 


wh«r« ® 1# Y t j, > 0, «• U. 

Let > 0 be such that. 


e t l l x J I 2 i v< Xl x N ,t) 

ut v o 2 v(x, (o) yo.o) 


If Vic^, 


i * l Vf) 

JeI 2l J 3 


V ^ 

Ox 2 


then V X 


v Y 


c(0, a - I y <t*) 2 ) 

JtI 2i 3 3 

V(*, V e> I - I x 1 1* | 

A " iC^ * 1 




V t > 0 


(2.5) 


( 2 . 6 ) 


2.4 Exponentially Stabilizing Control Laws 

In [1], various new exponentially stabilizing compensators were Introduced for both the set-point and 
tracking control problems. For the purposes of adaptive control, it is of interest to consider the subset of 
this class summarized In Table I. In addition, the well-known computed torque control has also been Included 
in Table I for comparison purposes. It is noted chat the desired potential field U*(Aq^) used in [1] has been 
chosen here simply as, 

U*(a<) 1 ) - Y Aqj 1 Kp^. <2. 7) 

so as not to obscure the presentation with additional obstacle avoidance objectives. Nevertheless, many of 
the adaptive control results presented herein are easily extended to the more general case. 

It is useful to observe that all Control Lavs 1-7 differ from the computed torque method in that the mass 
matrix M(q^) does not premult iply the position and velocity feedback galna K. and respectively. This 
property Is critical since It renders this entire class of control laws amenable to simple adaptation schemes 
(i.e., certainty equivalence adaptation) which can be shown to lead to desired asymptotic stability properties. 
The presence of the mass matrix premultiplier otherwise prevents simple cancellations In the Lyapunov function 
derivative, hindering most attempts to apply adaptive control directly to the nonlinear dynamic manipulator 
equations. A recent exception to this can be found in the work of Craig, Hsu and Sastry [9]. 

However, the resulting adaptation law requires that ths estimated mass matrix be invertible for all values of 
estimated parameters. This in turn requires on-line projections of parameter estimates Into prespecified 
bounded regions of parameter space where M(qj) is not only invertible, but where the true parameters are 
certain to lie. This approach not only requires tight bounds on parameter uncertainty, but involves a very 
difficult (al belt off-line) determination of the proper parameter projection domains. This problem is 
further exacerbated by the fact that the adaptation law is not parameterized by physical parameters and Is 
of the form where the transformation back to physical parameters is neither straightforward or unique. These 
problems are overcome in this paper by using the exponentially stabilizing control laws of Table I, which do 
not involve a premultiplying mass matrix on the feedback gains. 
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TABLE I 


COMPUTED TORQUE CONTROL LAW 


CONDITIONS FOR STABILITY + 

CROSS-REFERENCE 
TO (I]* AND (20| 

u * 

■ - H< qi )oyq, 

♦ tt y Aq 2 > + ktqj) ♦ ♦ C(qj.q 2 )qj 

None required 

(2.34) 

NEW 

EXPONENTIALLY 

STABLE CONTROL LAWS 




1. 

u - - K/q, - 

V 2 ♦ "‘V + M(a l )a 2d 

- $ J(qj.q 2 )q 2d ♦ 5 VV q 2d )q 2 

None required 

(4.2s) 

2. 

u • - *^<1, * 

* H(, i )q 2d * 

2 J( V , 2d >, 2 * 2 VV q 2 )q 2d 

None Required 

(4.2b) 

3. 

u • - K p A qi - 

Vq 2 ♦ k< V * M( V q 2d 

- Y J< V a 2d )q 2d * * VV q 2 ,q 2d 

None Required 

(4.2c) 

4. 

. - - K p iq 1 - 

V a 2 * k <V + H< V q 2d 

- } J( qi .q 2 )q 2 ♦ 1 VV q 2d )q 2 

None Required 

(4. 2d) 

5. 

* • - V q i • 

V q 2 * k(a l ) * M(, l >a 2d 

* C( V q 2d >q 2d 

WV * r 

(4.6) 

6. 

» - - V q ! * 

V 2 * k( V + H( V a 2d 

♦ C(q r q 2 )q 2 

a sin**v* Bufflcl * ntl y l * r R* 
w.r.t. Initial condition 

(4.7) 

7. 

u . - K p iq l - 

X y iq 2 ♦ k(q ld ) ♦ H(q u )q 

2d * C(q ld ,q 2d )q 2d 

o J (K ) sufficiently large 
■ in v 

w.r.t. Initiel condition 

(4.8) 


^General Aaiunpcloni lh 2 d ll. I I I I boundtd ; k y * * °* K p " * p T * 0 

*Ut UMAdj) - j Aq^ K p A qi in (1| 


In the nonadaptlve case, comparisons between chs ntv control laws of Tabls I and ths computsd tor qua 
msthod can be found in (!]. Nevertheless, a brief account Is in order here. In particular. Control Lews 1, 
2,3,4 are roughly "on par" with the computed torque method in the nonadaptlve case, guaranteeing exponential 
stability with no conditions on Kp or K v . Unlike the computed torque method, however, they are not in s form 
suitable for application of the recursive Newton-Euler computation technique. This presently appears to be 
their major disadvantage. In order to overcome this difficulty. Control Laws 5, 6 and 7 were developed 
in a form suitable for recursive Newton-Euler computation. Relative to the computed torque method. Control 
Law 5 utilizes the desired velocity signal q 2d in place of the measured velocity q2 In the nonlinear terms of 
the controller. This "cleans up" the feedback signal in the sense that nonidealities due to sensor dynamics 
and measurement noise in q 2 are avoided in the nonlinear feedback terms. Control Law 7 further replaces q^ 

In K, M and C by qj<j. This decouples the nonlinear terms from real-time measurements, which completely 
removes the requirement for on-line computation of nonlinear terms in the controller implementation. Control 
Law 6 is exactly the computed torque method without the premultiplying mass matrix term described earlier. 

The advantages of these controllers are off-set slightly by the conditions imposed on Kp and K v fo* guaranteeing 
asymptotic stability i.e., that K v be chosen sufficiently large for Control Laws 1, 2, 3, 4, 5, 6 and that both 
K v and K p be chosen sufficiently large for Control Law 7. It will be seen in the adaptive case that these 
requirements can be removed by adapting these feedback gains appropriately. 

The use of q2<j rather than q 2 in many of the new control laws offers additional advantages. In particular, 
in the set-point control application q2d"^2d“ 0 ' Henc *» there is considerable simplification in the control 
laws relative to the computed torque method, i.e., the nonllneer terms vanish from the control law. This 
simplification carries over directly to the adaptive case and provides substantial simplification in set-point 
control relative to the recent adaptive control laws of Slotlne and Li (11] and Paden(10J. 


3. A New Class of Asymptotically Stable Adaptive Control Laws 

All of the new exponentially stabilizing control laws suuaarlzed in Table I have the unique 
property that can be adapted in real-time so as to yield asymptotically stable adaptive control systems. 
Furthermore, the adaptation is done in a certainty equivalence fashion, I.e., by simply replacing 
unknown quantities in the control laws by their estimates - as generated by an appropriate parameter 
adaptation algorithm. In this section, asymptotic stability for ths various control lavs will be 
proved, and the proper mechanisms for parameter adaptation will be derived. 
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The simplicity la ttructurm of the adaptive control schemes presented horo la largely duo to a 
*1 inoar In tha parameters" formulation of tha problem. This particular paraaatarlsatlon la bacoalng 
Increasingly popular In racant lltaratura (cf., [9] (10] (16) [19D and will bo dlacuaaad In nora detail 
bo low. 

3.1 Llnaar In tha Paraaatara Formal at ion 

A uaaful paraaatarlsatlon of tha nonlinear dynamical aquation* arlaoa by noting tha following 
relatione (x, y and t arbitrary vectors), 

C(x,y,)y • H (x,y)e 
C C 

M(x)s • H M (x,s)e H 
k(x) * 

V x>y)y - HpCx.yJdjj 

where He, Hr, H* and Hj> are known matrix valued functions of x, y and s, and where « c , 0 m* ®k» ®D 
vectors of constant parameter* related directly to true physical paraaatara (massej, Inertias, link 
length*, center of gravltlea, etc.). It la emphasised that this paraaatarlsatlon does not contain 
any hidden "slowly varying" states fn tha parameter vector definition and does not require any 
linearisation of tha dynaalcal aquations of aotlon. 

3.2 Clobal Asymptotic Stability for Adaptation of Control Lavs 1, 2, 3, 4 


In this section, global asymptotic stability Is proved for adaptation of Control Laws 1, 2, 3 and 4. 

In order to avoid redundant analysis, tha details of tha proof will ba considered only for Control Lav 1, 
and the extension to tha other control laws will follow immediately by taking advantage of the unified 
treatment of these control laws given In (1]. 

3.2.1 Asymptotic Stability 

Consider Control Lav 1, 

u°- - K p 4vV<2 * ' I J( V q 2 )q 2d I W^d 1 ^ (31) 

Here, superscript la used to denote the ideal nonadaptlve control lav, l.e., the completely "tuned" 
control law which would be used if the parameters were known exactly. Using the linear in the parameters 
formulation discussed in Sec. 3.1 there exists a matrix H 1 (qj, q2, q2d» 92d> and * v<ccor °* parameters 
8 such Chat, 

u°- - K p A<l 1 - r v A<, 2 + H t 8 (3.2) 


where 


V * "<1x>i 2 d - 1 J(q l* q 2 ), '2d + 1 M D (<, l ,q 2d )<, 2 ° >3> 

Here, the parameters in ® are constant with time and are related directly to physical link and payload 
parameters . When these parameters are unknown, the parameter vector & is replaced by Its estimate 
8(t) In real-time to give the following adaptive control law, 

u • - K - K,Aq 2 * P- 4 > 

Subtracting (3.2) froa (3.4) and rearranging glv.a 

u • u° + Hj ( 8 - 0 ) - u° + Hj* ( 3 . 5 ) 

This is an Important relation alnce it shows that the adaptive control la equal to the nonadaptlve control 
plus an expression which is linear in the parameter error * * 0-8. 

The proof of stability then followa by choosing the following Lyapunov function, 

V - V° + J * T 1*4 r - r T > 0 ( 3 . 6 a) 

T 

where V° is the Lyapunov function for the nonadaptlve control law used In [1], and where ♦ T* Is a positive 
definite function in the parameter error > ♦ For completeness, V° Is rewritten here (cf., [1],(4.4) where 

U*(Aq,) - J Aq, T K^Aq,), 
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(3.6b) 


v ° ■ I S* + i \ T + c V "(VS 


Taking tha derivative of V along ays t am t raj act or la a and aubatltutlng control lav (3.3) given upon 
rearranging, 

V - V° ♦ (6q 2 ♦ ci^ l ) T Hj 6 + i T l> ( 

where V° la tha Lyapunov function darlvatlva for tha nonadaptlvo cane, and- where the additional tarva 
Involving 4 on tha right hand alda of (3.7) arlaa dlractly fron tha additional tarva Involving 4 In tha 
control lav (3.3) and tha Lyapunov function (3.6) raapactlvaly. 

Tha aacond and third tarva of (3.7) ara cancelled exactly by tha choice of adaptation lav. 


i - e - - r _1 B 1 T (6q 2 ♦ c4q x ) 


Tha .xprasalon (or tha raaalnlng t.ra V° 1. alaply takan (roa (1) aa, (aaa (1), (4.5) vhara v ft o >ln (Kp), 
also note that Control Lav 1 correaponda to caae (4.2b) for which a * 2 ) 

2 


« X IKJI 2 - « 2 l|4q 2 H 2 + T 21 IN 11 HMjH 2 


(3>il 

°2 *WV (3 - 1( 

Y 21 • f \ (31( 

h, “ aa* I lq,.l |h (3.1) 

2 ’2d 2d 1 

h, ■ aa* (J | |n (q )| |) (3.1: 

1 ’1 1 11 

w - aas ||M(q )|| (3 . 1J 

q l 2 

0 < c < 1 arbitrary 
o* arbitrary 

t* arbitrary 

Applying tha 8-ball arxuaant of lama 2.1 to (3.9) ualng the valuea of at, a?, and y-m given in 
(3.10), It followa that If, 


. 3 3 

a . (K ) > t n- o 
•in p 4 2 


, n, n, V -r 

WV * e <- ♦ r-f ♦ t *r> > 
0 1 


»:■ ‘illiojll 2 - * 2 IIa« 2 II 2 


for any X ^ and \ 2 auch that. 


X lt (°. c(a atn « p ) -|n 2 p 2 )) 1 

x 2 =»• WV- c( ^ : i + r <r> J > 

0 1 


V o • V | 


C 1 • 7 t 0 -in% + <V - t2 c 2< M >) 
, 2 . i a . £_) a 00 

o (M) ft Bin o (M(q. ) ) 

“ ln 1 
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Since p Is Arbitrary, It can bs chosen sufficiently small so that (3.12a) la satisfied. Vlth this 
choice of P 2 , the value of c in (3.12b) can be chosen sufficiently nail so that Inequality (3.12b) Is 
satisfied. Hence (3.13) follows. This Is essentially the sane proof of stability as In the nooadsptlve 
case (c.f., (1] Theorem 4-1) with the following exceptions, 

1) The value of V Q In (3.12b) and (3.14b) now Includes the Initial parameter error j 4 T (0)r+(0). 

2) The value of c Is now required for Implementation of the parameter adaptation law (3.S). 

3) ♦ In (3.13) Is now only negative senldeflnlte In the state since the full state vector in the 
adaptive case Is augmented by 4. 

It Is noted that proporty 3 destroys the simple exponential stability ariument used earlier In the 
nonadaptlve case (cf., (lj, Theorem 4-1) to Insure asymptotic convergence ofllAqJI and I|AqJ|. In addition 
since tbs error eystem In (Aqt, Aqj»4) Is oonautooooous (and In general, nonperiodic), standard Invariance 
principles are not applicable. Alternatively, we make use of a leaa due originally to 
larbalat 9 quoted without proof from Fopov [14] (pg. 211^. 

Immt 3-1 (Barbalat) 

If V Is a real function of the real variable t, defined and uniformly continuous for t > 0 and if the 
limit of the integral ^ 

11m f W(t*)dt' , 

t-M * O 

exists and la a finite number, then 

11a V(t) - 0 . ■ 

t— 

For our purposes let, 

V(t) $ X l ||4q 1 (C)|| 2 + X 2 ||4q 2 (t)|| 2 
so that 

V<-w (3-^7) 

Integrating both sides of (3.17) from 0 to t, yields upon rearranging, 

, t 

«Jt' < » o • V(t) (3.18) 

o 

Since V 0 Is bounded, and V(t) Is nonincressing and bounded below, It follows that 

11a f W dt' < • 
t— o 

Also, since W is bounded, W(t) is uniformly continuous. Hence, application of Barbalat’s lam a gives, 

11a W - 0 (3*19) 

or equivalently | | 1 1 -^0 and | | 1 I 

This completes the proof of asymptotic stability. The proof, however, is not a global one due to 
property 2, l.e., the value of c which was not required in the nonadaptlve case now appears in the 
parameter adaptation law (3.8). Hence, one Is cOMltted to choosing a particular value of c in the 
adaptive implementation. Of course, c can always be chosen sufficiently small to satisfy the requirement, 
however, the position tracking performance determined by the magnitude oZ lj in (3.14a) must be compromised 
as a result. Bence In practice, the initial choice of c can be made using whatever bounds on ni» 02« u» 
o(M) and V 0 are available a-prlorl, and the value of c can be Improved (increased) on-line ss more information 
becomes available. It Is noted that (3.16a) and (3.16b) impose additional constraints on how large c can 
become, since it is required chat > 0 and > 0 for * positive definite V (these conditions can be shown 
sufficient) . 

The asymptotic stability proof presented above for adaptation of Control Law 1, is easily extended to 
adapt at ion of Control Lave 2, 3 and 4, since the corresponding nonadaptlve Lyapunov function derivatives 
for these control levs are of exactly the same form aa V® in (3.25) (see [1], Theorem 4-1 for details). 

For convenience, all asymptotically stable adaptive control. lavs discussed thus far, and their appropriate 
parameter adaptation lavs are summarized in Table II, corresponding to cases l.s, 2.s, 3. a, and 4. a, 
respectively. 

An alternative to choosing c sufficiently small in the above asymptotic stability argument is to 
choose K v sufficiently lsrge. In this esse, the condition on c above can be removed completely by adapting 
ILy on-line. This modification insures global asymptotic stability of the adaptive control system (l.e., choice 
of c Independent of the initial condition V Q ) and is discussed in more detail below. 
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3.2.2 Global Asymptotic Stability-Adapting Ky 

Sinca the velocity gain Ky enters linearly in the control law, it can be adapted as if it were an unknown 
parameter using the same formulation of Sec. 3.2.1. It will be shown that this approach removes the dependence 
of the choice of c on the initial condition V 0 and this completes the proof of global asymptotic stability 
for the adaptive case. 

Consider Control Lav 1 written in adaptive form, where both 0 and are adapted in real time i.e. f 

u - - y qi - v , 2 + V (3.20) 

litre, |L it a time-varying quantity which remains to be specified, and is as defined earlier in (3.3). 

The nonadeptive control lew u° in (3.1) is subtracted from (3.20) to give the following expression, 

u - u° - « T 4 <i 2 + Uji <3- 21) 

whar. AKy • ty ~ Ky .ad ♦ “ 

The Lyapunov function for the stability analysis is given as 

* - v° ♦ 7 * r r* * y « t*{a k v t « w ) . «>o, r-r T >o 0 . 22 ) 

where a new term has been added relative to (3.6a), quadratic in the error Taking the derivative 

of V along system trajectories end substituting control law (3.21) gives upon rearranging 

v - V° ♦ (Aq 2 ♦ cAq 2 ) T Hj# + pT* 

+ TR{J«Af v T - Aq 2 (A <) 2 + cAq^jA*,) (3.23) 


The latter terms are cancelled exactly by the choice of parameter adaptation lews, 

♦ - « - - r" 1 a x T (4q 2 + cAq : ) 

f ( 4 q 2 * c& tl^<*2 T 


(3.24.) 


(3.24b) 


The choice leaves V exactly of the form (3.9) i.e., applying the 0-Ball Leona 2.1, 
V-V°< - UlAq.ll 2 - X-l |Aq_; I 2 


* I n 2 ° 2 j < 

0 1 

In (3.26) and (3.27), all quantities are defined exactly as in (3.12a) and (3.12b) respectively, except 
for V Q which is presently the initial value of V in (3.22). Furthermore, the values of and are 
once again given as 


*1 * I 1 °.ln (K p + C V * * 2 c 2 < m » 

e, - i u - £_ )« (m) 


An important observation is that. 


V o - <Xl|!C v in , IUC V 


CXllKjl) . UK 


Hence, for any choice of 6>0, c>0 and K_ ) - Kp^ > 0, there exist values of 0 , i^, and Ky*K v '^ > 0 (with °nin^v^ 
sufficiently large) such that inequalities (3.26) and (3.27) are satisfied, and Cj^O, $ 2 *® * Q 0*28) end 
(3.29), respectively. Global asymptotic stability of this adaptive control scheme then follows immediately 
by application of Barbalat's Learns to the Lyapunov function derivative (3.25), as was done earlier in 
equations (3.17) through (3.19). 
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Tba global asymptotic stability of adaptive controllers bawd on Control Lava 2, 3 and 4 (whara Ky la 
adaptad on-line) follow from an idantical argument, slnca 4° corresponding to tba nonadaptlva Lyapunov 
function darivatlvaa for thaaa control lava ara of exactly tba aa me fora aa $° In this ana ly a la (saa (1], 
Theorem 4*1 for details}. for convenience, thaaa adaptlva control lava involving adaptation of K* ara 
summarised in tabla II, corraapondlng to caaaa l.b, 2.b, 3.b, and 4.b, raapactlvaly* 


3.3 Global Asymptotic Stability for Adaptation of Control Lava 3, 6 and 7 

Global asymptotic atabillty for adaptation of Control Lava 3, 6 and 7 can ba proved ualng exactly tba 
ttae technlquea aa applied in Sac. 3.2. The only dlffaranca Ilea in alight varlatlone in the nonadaptlva 
Lyapunov function derivative V° which ariaaa In each adaptlva control analyela 

Due to apace limitations, thaaa proof a have bean oalttad, but the raaulta ara summarized In Tabla II 
corraapondlng to caaaa 3. a, 5.b, 6. a, 6.b, and 7. a, 7.b, raapactlvaly. Dataila can ba found in [21] , to 
which the aquation numbers In Tabla II ara referenced. 


4. Adaptlva Computed Torque Method 

It vaa Motioned earlier that In the computed torque Mthod (l.a. v control law 0) the praaanca of the M(q^) 
tarn premultiplying the Kp and Ky galna complicates the Lyapunov analyala and hlndara moat alapla attempts to 
Mka It adaptive. Nevertheless, the computed torque controller la a well-known control lav In the literature 
and la widely applied In practice. Hence, It la uaeful to Invcatigate under what condltiona It can be made 
adaptive, and to what extent adaptive atabillty can be guaranteed. For thle purpoae, we conelder a apeclal 
caae of the computed torque control lav which has acalar galna kp and ky, l.e., 

u° - - »(q 1 )(k p Aq 1 + k/q 2 ) + k(q x ) + Jtfq^q^ + c <V q 2 )q 2 


This la written In adaptive form aa, 
u * u° + Hg4 

where ♦ • 8 - 0 # and the linear in the parametera part has been chosen aa, 

H g 0 - - M(q 1 )(k p Ac ll *fk v Aq 2 ) + kCq^ + MCq^q^ + 


(4.2) 


(4.3) 


Let a Lyapunov function be defined as, 

v - v° + | ♦ T r$ , r - r T > o (4.4.) 

where 

v - \ 4q 2 T M( qi )aq 2 + | (k+ck^Aq^ M( qi )aq x + cAqJ M(q x )Aq 2 (4.4b) 


Then, the derivative of (4.4) along system trajectories Induced by control (4.2) Is given by 


V - Aq 2 T [-k v M(q 1 )Aq 2 + j ^(qj.Aq^qj) 


+ cAq x T l-k p M(q 1 )Aq 1 + M D (q 1> Aq 2 )q 2 + M^q^Aq^qjl 

(4.5) 

♦ (Aq 2 + cAq x ) T H g * + » T I> 

(4.5) 

Let, 


♦ - 9 - - r" 1 H g T (Aq 2 +cAq 1 ) 

(4.6) 

Then, 


V < - o 1 ||Aq 1 || 2 - “ 2 1 1 Aq 2 1 I 2 + Y 12 ||Aq 2 || | f Aq^ | | 2 


+ ^ 21 l|4q i l| + Y 22 l|Aq 2 ||)||Aq 2 || 2 

(4.7) 

where 


o - c(k o(M) - n,(l + o 2 )) 

1 p - 2 

(4.8) 

“2 ' V (M) " I n i | 
0 

(4.9) 

Y 12 * Cn l 

(4.10) 

v 2l ' cn l : t 22 ' I n l 

(4.11) 


239 


Applying the a -Ball Lemma, It follows that, 


V i- A 1 l|Aq 1 l| 2 - X 2 |U, 2 || 2 

(4.12) 

V ± 

ck o(M) » n,(l + p 2 )) + c n.Cr 2 ) 2 

p - 2 2 

(4.13) 


(4.14) 

c i “ i tk P * C K - cl2 )2< M >» 

(4.13) 

C 2 - | U - Sj)o(M) ; l 2 » 0 arbitrary 

(4.16) 


It Is noted that for any c > 0, both kp and ky can always ba chosen sufficiently large so that £i>0 and 
C 2 >0 (for appropriate choice of t 2 >0 in (4.15), (4.16)), and inequalities (4.13) and (4.14) are satisfied. 

Hence, the adaptive computed torque control law given by (4.2), (4.3) with parameter adaptation (4.6) is 
asymptotically stable when kp and ky are chosen sufficiently large. 

Since kp and k v must be chosen sufficiently large with respect to the Initial condition V 0 (c.f., (4.13), 
(4.14)) this proof of asymptotic stability Is not global (i.e., for fixed k p and ky there will always exist some 
V 0 such that *i and/or *2 ar « not positive). For this particular algorithm, it is presently not clear how to 
adapt k p and ky to Insure global asymptotic stability since the control u in (4.1) is not linear in the 
parameters (6, k p , ky) . 


5* Summary and Remarks 

The adaptive control lavs derived herein, along with the sufficient conditions for stability and appropriate 
parameter adaptation lavs are summarized In Table II. Several remarks are In order at this point in the 
discussion. 

Remark 5-1 All adaptive control laws in this paper were derived for the general tracking control lav. However, 
significant simplification occurs in many of these designs for the special case of set-point control (i.e., 
^2d“^2d“°) * 

Remark 5-2 The adaptive robustness issue remains open. Certainly for parameter adaptation laws of the form 
given in Table II , there will be sensitivities to noise disturbances and unmodelled dynamics directly analogous 
to chose which arise in the linear adaptive control case. It presently appears that many of the robustness 
techniques developed in the linear adaptive control literature will carry over to the nonlinear adaptive control 
application. This conjecture, however, remains to be investigated. 

Remark 5-3 In the nonadaptive case, many of the control laws in Table II are In a form appropriate for 
application of the recursive Nevton-Euler computational algorithm. However, the Newton-Euler algorithm requires 
knowledge of all physical parameters-more parameters than^are actually needed to control the system adaptively 
and more than are actually adapted on-line In the vector 0 of Table II. Hence, the transformation from 0 back 
to physical parameters is required in order to salvage use of the Hevton-Euler algorithm in the adaptive case. 
However, the transformation is generally nonlinear and will not lead to a unique solution unless further 
constraints are imposed. One typical set of constraints arises when only the payload mass is unknown. In 
the more general adaptive case, it is useful Co note Chat all linear in the parameters expressions can be 
implemented directly, since representations of the form H8 are assumed to be available in symbolic form. 

Remark 5-4 The control laws of Table I were derived In [1] for the general desired potential energy function. 

This feature was dropped in the adaptive case in order to simplify the analysis. However, it appears that 
the adaptive control laws developed herein can be extended to the more general esse and this line of research 
presently under Investigation. 

Remark 5-5 A brief comparison with the recent results Psden (10] and Slotlne and Li (11] is useful. In (10] (11], 
adaptive control laws are derived by choosing u to cancel various terms in the Lyapunov function derivative, 
rather than overbounding them (via Lemma 2.1) as was done here. This approach has the advantage of providing 
global asymptotic convergence without adapting gains and Kp. The control laws, however, are by necessity more 
complex than those designs considered here, and do not simplify in the set-point control case. 

6. Conclusions 

A new class of asymptotically stable adaptive control laws is defined by adapting the control laws of (1] 
in a certainty equivalence fashion. These algorithms are proved to be asymptotically stable without approximations, 
linearizations or ad-hoc assumptions concerning the nonlinear manipulator dynamics. Furthermore, the asymptotic 
convergence properties can be made global by appropriate adaptation of feedback gains. On-going research 
efforts are directed at adaptive robustness, computation, and obstacle avoidance problems. 
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TABLE II SUMMARY OF ASYMPTOTICALLY STABLE ADAPTIVE CONTROL LAWS 




Moo. required Global A. W to etc Stability 


R 



Kacjrsiv* Mavton Eultr applicable In nonadapti*a cm 


Acknowledgements 


The authors ir« Indebted to Dr. K. Krauts for many helpful technical dlecusalons. This work was done 
at the Jet Propulsion Laboratory, California Institute of Technology under contract with the National 
Aeronautics and Space Administration. 

Inferences 


[1] J.T. Wen and D.S. Bayard, "Simple Robust Control Lavs for Robotic Manipulators - Part I: Nonadaptlve 

Casa," this Workshop. 

[2] S. Dubovsky and D.T. Desforges, "The Application of Model-Referenced Adaptive Control to Robotic 

Manipulators Trans. ASMS J. Dynamic Systems, Measurement and Control, Vol. 101, 1979, pp. 193-200. 

[3] A.J. Kolvo and T.-H. Guo, "Adaptive Linear Controller for Robotic Manipulators," IEEE Trans. Auto. 

Contr., Vol. AC-28, No. 2, 1983, pp. 162-171. 

[4] C.S.G. Lee and M.J. Chung, "An Adaptive Control Strategy for Mechanical Manipulators," IEEE Trans. 

Auto. Contr., Vol. AC-29, No. 9, 1984, pp. 837-840. 

[5] M. Tomlxuka and R. Horowitz, "Model Reference Adaptive Control of Mechanical Manipulators," IFAC 

Workshop on Adaptive Systems In Control and Signal Processing, San Francisco, CA, 1983. 

[6] M. Takegakl and S. Arlmoto, "An Adaptive Trajectory Control of Manipulators," Int. J. Contr., 

Vol. 34, No. 2, 1981, pp. 219-230. 

[7] K.Y. Lin and M. Eslaml, "Adaptive Controller Designs for P.obotlc Manipulator Systems Using Lyapunov 

Direct Method," IEEE Trans. Auto. Contr., Vol. AC-30, No. 12, 1985, pp. 1229-1233. 

[8J T.C. Hsla, "Adaptive Control of Robot Manipulators - A Review," IEEE Int. Conf. Robotics and 
Automation, San Francisco, CA, 1986. 


(*] 

[ 10 ] 

fllj 

( 12 ] 

(131 

(14] 

(15] 

(16] 


J. Craig, P. Hsu, and S. Sastry, "Adaptive Control of Mechanical Manipulators," IEEE Conf. Robotics 
and Automation, San Francisco, CA, April 1986. 

B. Paden, "PD+- Robot Controllers: Tracking and Adaptive Control," Private Communication. 

J. -J.E. Slotlne and W. Li, "On the Adaptive Control of Robot Manipulators," ASME Winter Meeting, 

Anaheim, CA, 1986. 

I.D. Landau, "A Survey of Model Reference Adaptive Technlquee - Theory and Applications," Automatics, 
Vol. 11, 1974, pp. 353-379. 

K. S. Narendra and R.V. Monopoll, Eds., Applications of Adaptive Control , Academic Press, New York, 1980. 

V.M. Popov, Hyperatablllty of Control Systems , Springer-Verlag, New York, 1973. 

P.A. Ioannou and P.V. Kokotovlc, "Robust Redesign of Adaptive Control," IEEE Trans. Auto. Contr., 

Vol. AC-29, No. 3, 1984, pp. 202-211. 

P.A. Ioannou, "Robust Adaptive Controller with Zero Residual Errors," IEEE Trans. Auto. Contr., 

Vol. AC-31, No. 8, 1986, pp. 773-776. 


[17] D.S. Bayard, C.H.C. Ih and S.J. Wang, "Adaptive Control for Flexible Space Structures with Measurement 
Noise," submitted to American Control Conference, Minneapolis, Minnesota, June 10-12, 1987. 


[18] C.C. Atkinson, C.G. An and J.M. Hollerbach, "Estimation of Initial Parameters of Manipulator Loads and 

Links," International Symposium on Robotics Research, 1985. 

[19] P. Khosia and T. Kanade, "Parameter Identification of. JtohtpJL Dynamics." IEEE Conf. Decision and Control, 

Fort Lauderdale, Florida, 1985. ' - - 

[20] J.T. Wen and D.S. Bayard, Robust Control for Robotic Manipulators, Part I: Noil-Adaptive Casa,(JFL 
Internal Document, Engineering Memorandum, EM 347-87-203). Jet Propulsion Laboratory, 

Pasadena, California, 1987. 

[21] D.S. Bayard and J.T. Wen, Robust Control for Robotic Manipulators, Part II: Adaptive Case, (JPL Internal 
Document, Engineering Memorandum, EM J4 7-87-204). Jet Propulsion Laboratory, Pasadena, California, L987. 


243 



